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Abstract 

The fidelity and entropy in an easy-axis antiferromagnetic Heisenberg spin-1 chain are studied numer- 
ically. By using the method of density-matrix renormalization-group, the effects of anisotropy on fidelity 
and entanglement entropy are investigated. Their relations with quantum phase transition are analyzed. 
It is found that the quantum phase transition from Haldane spin liquid to Neel spin solid can be well 
characterized by the fidelity. The phase transition can be hardly detected by the entropy. But it can be 
successfully detected by the first deviation of the entropy. 
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I. INTRODUCTION 



In condense matter physics, quantum phase transitions imply fluctuations, which happened at 
the zero temperature |T]. When a controlling parameter changes across critical point, some prop- 
erties of the many-body system will change dramatically. Many results show that entanglement 
existed naturally in the spin chain when the temperature is at zero. The quantum entanglement 
of a many-body system has been paid much attention since the entanglement is considered as 
the heart in quantum information and computation [21 As the bipartite entanglement mea- 
surement in a pure state, the von Neumann entropy in the antiferromagnetic anisotropic spin 
chain [5] and isotropic spin chain |6] are investigated respectively. By using the cross fields of 
quantum many-body theory and quantum-information theory, von Neumann entropy is applied to 
detect quantum critical behaviors [714TT]. A typical example is that Osborne solved exactly one- 
dimensional infinite-lattice transverse-field Ising model to obtain entropy by the Jordan- Wigner 
transform. The entropy predicts the quantum phase transition successfully [S] . Moreover, another 
concept from quantum information science, the ground state fidelity has been used to qualify quan- 
tum phase transitions in the last few years [ISHSlj. Because the fidelity is a measure of similarity 
between states, the fidelity should drop abruptly at critical points as a consequence of the dramatic 
changes in the structure of the ground states, regardless of what type of internal order is present 
in quantum many-body states. This result is the orthogonality of different ground states due to 
state distinguishability. Many results shown that the fidelity and the entanglement entropy have 
similar predictive power for identifying quantum phase transitions in the most systems [231 El]- 
However, the ground state fidelity is a model-dependent indicator for quantum phase transitions. 
It cannot be used to characterize the Berezinskii-Kosterlitz-Thouless (BKT) transition [251 126] . 
The BKT-like transition occurs at Heisenebeg isotropic spin chain with next-nearest-neighbor and 
in antiferromagnetic anisotropic Heisenberg model with A = 1 [23]. Entanglement, fidelity and 
their relations with quantum phase transition in high spin chain, just like spin-1 chain need to be 
further investigated. 

In this paper, the fidelity and entanglement in the easy-axis spin-1 chain are numerically inves- 
tigated by using the density-matrix renormalization-group (DMRG) technique. In section II, the 
Hamiltonian of easy-axis spin-1 chain and its critical property are presented. In section III, the 
effects of anisotropic interaction on fidelity is investigated and its relation with quantum phase 
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transition is analyzed. The effects of anisotropic interaction on entanglement entropy is calculated 
and its relation with quantum phase transition is analyzed in section IV. At last, a discussion 
concludes the paper. 

II. HAMILTONIAN AND ITS CRITICAL PROPERTY 

It is known that there is a gap in the spectrum of isotropic Heisenberg spin-1 chain in the 
thermodynamic limit. The gap between the ground state and the first excited state energy is 
usually called Haldane gap [27]. The Hamiltonian of an anisotropic Heisenberg antiferromagnetic 
spin-1 chain of sites can be given by 

N-1 

H = jY,iS!S!^, + S^Sl, + AS^S^,), (1) 

i=l 

where S°'{a = x, y, z) are spin operators on the z-th site and is the length of the spin chain. The 
parameter J > denotes the antiferromagnetic coupling and J = 1 is considered in the paper. The 
parameter A is anisotropic interaction, < A < 1 is the easy-plane anisotropy and A > 1 is the 
easy-axis anisotropy. It is predicted that a novel phase Aci G [0 ~ 0.2] < A < Ac2 = 1-17 appears 
between the XY phase and the Neel phase [28] - [30iJ . Many results suggested that a BKT transition 
occurs at Aci-, while that the transition at belongs to the 2D Ising universality class [311 ^ 
[51] . Since the fidelity cannot detect a BKT-like phase transition [2S1 ISE] , we only concentrate on 
quantum phase transition which happens at easy-axis anisotropy. It is found that the Hamiltonian 
in Eq. (1) can be transformed conveniently to give 

N-1 

H = jY, [P(^r^l+i + S^Sl,) + SlSt^^ , (2) 

i=l 

where p = 1/A. The critical point of transition from Haldane spin liquid to Neel spin solid is 
around pc — 0.85. It is noted that the model of Eq. (2) is invariant for p — )■ —p, so that the results 
in the paper can be equally applied to a system with a ferromagnetic coupling on the xy plane. 

III. FIDELITY SUSCEPTIBILITY 



The ground state fidelity can be applied to detect the existence of the quantum phase transitions. 
A general Hamiltonian of quantum many-body system can be written as H{X) = Hq + XHj where 
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Hi is the driving Hamiltonian and A denotes its strength. If p(A) represents a state of the system, 
the ground state fidehty between p(A) and p(A + 5) can be defined as 

F(A, 5) = Tr[^pV2(A)p(A + 5)pi/2(A)]. (3) 

If the state can be written as p = \ip) theEq. (3) can be rewritten as F(A, 5) = \{ip{X)\ip{X+6))\. 
Because -F(A, 6) reaches its maximum value Fmax = 1 for 5 = 0, on expanding the fidelity in powers 
of (5, the first derivative ^^^a'f~°^ = 0. By using the property, the fidelity can written by 



F(>.J)^l + m^h.,S\ (4) 



Therefore, the average fidelity susceptibility ^(A, 6) can be given by [HI [19 



Q(x r 2[1-F(A,^)] 

It is well known that it is hard to calculate the ground state fidelity because of the lack of 
knowledge of the ground state function. For models that are not exactly solvable, most of re- 
searchers resort to exact diagonalization to obtain the ground state for small size. This method 
cannot precisely quantify the quantum phase transition because the size of the system is too small. 
Recently, the method of density-matrix renormalization-group (DMRG) [351 136] can be applied 
to obtain the ground state of the model. Moreover, the technique of calculating the overlap of 
two different ground states by DMRG has been used [371 EH]- For high precision, the method 
in [Sn] is used to calculate the ground state fidelity susceptibility. We calculate up to 80 when 
S = 0.001 [23j. The total number of density matrix eigenstates held in system block is m = 70 in 
the basis truncation procedure. The Matlab codes of finite size density-matrix renormalization- 
group have double precision. They are performed with three sweep in private computer, and the 
truncation error is smaller than 10"^''. 

The ground state fidelity susceptibility 5* is plotted as a function of anisotropic parameter p for 
different sizes in Fig. 1 in dimensionless units. It is shown that there is a peak in S. The peak 
of S increases when the size increases. The location of the peak decreases slightly to small value of 
p as N increases. The ground state fidelity F is plotted as a function of p for different sizes in the 
inset of Fig. 1 in dimensionless units. It is seen that there is a sharp valley corresponding to 
the peak of S. The critical point of fidelity susceptibility 5* and fidelity F is about pc = 0.85. As we 
known, the fidelity measures the similarity between two states, while quantum phase transitions are 
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intuitively accompanied by an abrupt change in the structure of the ground state wave-function. 
This primary observation motivates researchers to use fidehty to predict quantum phase transitions. 
Up to now, the ground state fidehty and fidehty susceptibihty have been apphed in various many- 
body systems to detect quantum phase transitions successfuUy [T2V[2^ . It confirms further that 
the Haldane spin hquid— Neel spin sohd transition occurs at the point in our system. 

IV. ENTANGLEMENT ENTROPY 

For comparison, the ground state entanglement entropy is also used to detect the quantum 
phase transition. For density matrix p of any pure state, the entanglement entropy Eab between 
subsystem A and B can be defined as 

Eab = -Tr{pA log2 Pa) = -Tr{pB logs Pb), (6) 

where pA{B)is the reduced density matrix obtained from p by taking the partial trace over the state 
space of subsystem B{A). 

In order to avoid boundary effects, the entanglement entropy of two neighboring central sites of 
large system needs to be calculated by using the method of DMRG. The entropy of two neighboring 
central sites labeled E]\f/2, N/2+1 is plotted as a function of anisotropy p with sizes of = 40, 60, 80 
in Fig. 2 in dimensionless units. The entanglement entropy of the two neighboring central sites 
increases monotonously when the anisotropy p increases. It is known that the peak or discontinuity 
in entanglement entropy indicates the quantum phase transition. The peak does not occur at the 
quantum critical point. It seems that this is due to the monogamy property [9j. We also calculate 
the entanglement entropy between the right-hand N/2 contiguous qubits and the left-hand N/2 
contiguous qubits. The result is similar with that shown in Fig. 2. For avoiding repetition, the 
result is not plotted again. This means that the entropy cannot be used to predict quantum phase 
transition here. However, the critical properties can be captured by the derivatives of the entropy 
as a function of the anisotropy p [9l HHl HI] . 

It is noted that the transition from Haldane spin liquid to Neel spin solid is continuous and 
belongs to a second order quantum phase transition. It has been shown that the density matrix may 
be continuous but its first-order derivative is singular at a second order quantum phase transition 
point [12]. Consequently, we calculate the first-order derivative of entropy dE]\r/2^N/2+i/dp at the 
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ground state with 6 = 0.001. The first derivative of the central two sites entanglement entropy 
dEN/2,]Sf/2+i/dp is plotted as a function of anisotropy p for different sizes in Fig. 3. It is seen that 
there is a peak in dEN/2,N/2+i/dp. The location of the peak in dE]sr/2,N/2+i/dp moves slightly to 
small value of p as increases. The critical value is about Pc — 0.85. It represents indeed the 
quantum phase transition from Haldane spin liquid to Neel spin solid. 

V. DISCUSSION 

In the paper, the fidelity susceptibility and entropy in an easy-axis antiferromagnetic Heisenberg 
spin-1 chain are studied. By using the density-matrix renormalization-group for the model, the 
effect of anisotropic interaction on fidelity susceptibility and entropy of large size is presented. 
Their relations with quantum phase transition are investigated. It is shown that the quantum 
phase transition from Haldane spin liquid to Neel spin solid is clearly marked by the peak (or 
valley) of the fidelity susceptibility (or fidelity). However, the entanglement entropy cannot have 
similar predictive power for revealing quantum phase transition in the system due to the monogamy 
property, while the first-order derivation of entropy can be used to successfully detect the quantum 
phase transition. 
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FIG. 1: The fidelity susceptibility S of the ground state is plotted as a function of anisotropy p for different 
size N. The fidelity F is plotted in the inset. The symbols are for N = 40(B), iV = 60(»), iV = 80( 0- 
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FIG. 2: The two neighboring central sites entanglement entropy N/2+1 is plotted as a function of 
anisotropy p for different size N. 
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FIG. 3: The first derivative of tlie two neigliboring central sites entanglement entropy dEj^ j2^N/2+i/ is 
plotted as a function of anisotropy p for different size N . 
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